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Matrix models of Yang-Mills type lead to an emergent gravity theory, which may not require 
fine-tuning of a cosmological constant. We find cosmological solutions of Friedmann- Robertson- 
Walker type. They generically have a big bounce, and an early inflation-like phase with graceful 
exit. The mechanism is purely geometrical, no ad-hoc scalar fields are introduced. The solutions 
are stabilized through vacuum fluctuations and are thus compatible with quantum mechanics. This 
leads to a Milne-like universe after inflation, which appears to be in remarkably good agreement 
with observation and may provide an alternative to standard cosmology. 
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Quantum field theory and general relativity provide 
the basis of our present understanding of fundamental 
forces and matter. However, there is up to now no sat- 
isfactory way to reconcile them in a consistent quantum 
theory. General arguments based on quantum mechanics 
and general relativity suggest a "foam-like" or quantum 
structure at the Planck scale 10^^'^ cm. This problem 
has its most dramatic manifestation in the cosmological 
constant problem: the small but non-vanishing cosmo- 
logical constant in the currently accepted ACDM model 
is in striking contradiction with quantum mechanical ex- 
pectations, which are off by at least 60 orders of magni- 
tude. No satisfactory solution of this problem within the 
conventional frameworks has been found. 

We point out in this letter that emergent gravity on 
non-commutative (NC) spaces may provide a resolution 
of these problems. The starting point are matrix models 
of Yang-Mills type 



5- 



YM 



-TT[X%X'][X\X'>]rjaa'Vbb' 



(1) 



supplemented by suitable ferniionic terms. Such mod- 
els have been discussed in the context of NC gauge 
theory and string theory. Here a = 1,...,D are 
(infinite-dimensional) hermitian matrices, and rjab = 
diag(l, 1, ...,±1) in the Euclidean resp. Minkowski case. 
The basic hypothesis of our approach is that space-time is 
realized as 3-l-l-dimensional NC "brane" solution of such 
a model. The effective geometry and gravity on such a 
brane was clarified recently P, [H, : the effective metric 
is not fundamental but depends on the embedding and 
the Poisson structure, reminiscent of the open string met- 
ric . This leads to an emergent gravity closely related 
to NC gauge theory, as anticipated in [Sj, ;6;] . Among these 
models, the IKKT model is singled out by maximal 
supersymmetry, required by consistency at the quantum 
level. This implies in particular D ~ 10, which we as- 
sume from now on. This model is a candidate for a theory 
of all fundamental interactions and matter. 
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The model admits 4-dimensional noncommutative 
spaces M.e C. as solution, interpreted as space-time 
embedded in 10 dimensions. This can be seen by splitting 
the matrices as 

X'' = {X^,<i>'), M = l,-,4, z = l,...,6 (2) 

where the "scalar fields" = (jf{X^) are assumed to 
be functions of X^. We only consider the semi-classical 
limit of such a quantum space, indicated by ^. Then 
Xt^ ~ is interpreted as quantization of coordinate 
functions on A4, define the embedding of a 4- 

dimensional subnianifold A4 C M.^°, and 



[Xt',X'']r^i0'"'{x), ^l,Iy=l,. 



(3) 



can be interpreted as Poisson structure on A4 . The phys- 
ical fields arise from fluctuations in the matrix model 
around such a background. Therefore they live only on 
the brane Ai, and there is no higher-dimensional "bulk" 
which could carry any propagating degrees of freedom, 
unlike in string theory or standard braneworld-scenarios. 
As shown in @,IH , the effective metric for all scalar, gauge 
and fermionic fields propagating on A4 is given by 



det 91, 



where 



(4) 
(5) 

(6) 

is the induced metric on M. G^i, is dynamical, depend- 
ing on the embedding fields 0* and 9^'^ . Standard embed- 
ding theorems imply that G^j/ can describe in principle 
the most general metric in 4 dimensions. Therefore the 
matrix model defines a theory of space-time and gravity 
coupled to gauge fields and matter. 

The matrix e.o.m. [X", [X'' ]]riaa' = can be writ- 
ten in a covariant manner as [3| 

(7) 



Ag0' = Agx'" = 0. 
This implies the covariant equation for 9^'^ 
G'"'V^(e"0-;) - e-'^G^.r^a^ 
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(8) 
(9) 
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which relates 9^'^{x) with the metric G^^ . Here V de- 
notes the Levi-Civita connection with respect to G^^. 
In principle, 9^'^ ^ breaks (local) Lorentz invariance. 
However, 6^^'^ does not enter explicitly the effective ac- 
tion to leading order, and its presence through higher- 
order terms may be below experimental limits if the 
scale of noncommutativity K^^c is high enough. In fact, 
this spontaneous breaking of Lorentz invariance leads to 
massless gravitons as discussed below. Note also that the 
matrix model defines preferred coordinates , which are 
not observable and not in conflict with observation. 

Equations ([7|) imply that the embedding M C MP is 
harmonic w.r.t. G^jy. A particularly interesting case is 
given by geometries with 

G^i/ = 9tiu- (10) 

It is not hard to see Q that this holds if and only if 
•q = e'^ , which (in the Euclidean case) is equivalent to 
-kO ~ ±9 where ★ is the Hodge star. Then ^ simplifies 
as 

= 0- (11) 

These are formally the free Maxwell equations in the 
background geometry Gfj,u- In particular, they have 
propagating massless solutions with 2 physical helicities. 
However, there are no charged fields under this "would- 
be" U{1) gauge field in the matrix model. Rather, these 
modes turn into the 2 physical degrees of freedom of 
gravitons. Writing — 9^^ + F^i, on a flat (Moyal- 
Weyl) background, the metric fluctuations are 

h — —G ,9'^'PF —G ,9^''PF ~G F 9P^ II 

which are nontrivial and Ricci-flat, i?py[G + /i] =0 0, @. 

Quantization and induced gravity: The quantiza- 
tion of the matrix model llj is defined by 

Z = y" dX"e~^^^'[^] (12) 

(omitting fermions for simplicity) . Since one cannot sim- 
ply add an explicit Einstein-Hilbert term, the model is 
highly predictive. The cosmological solutions given be- 
low add to the evidence that it may provide a (near-?) 
realistic theory of gravity, with great advantages for the 
cosmological constant problem. 

Consider a perturbative quantization of the matrix 
model around a given background as discussed above. 
Since all fields couple to G^i^, standard considerations 
imply that in the classical- geometric limit, the effective 
action at one-loop can be obtained from Seeley-de Witt 
coefficients, so that 

ri-ioop = J^J d^^VWl (ciAt + C4i?[G] + 0(ln A)) . 

(13) 

The coefficients Ci as well as the effective cutoffs Aj de- 
pend on the detailed field content of the model, cf. 
This is essentially the mechanism of induced gravity. 



Now consider the equations of motion for the geome- 
try, taking into account the quantum contribution (|13p . 
Remarkably, equation ([5]) for 9'^" is unchanged: it is a 
direct consequence of Noethcr's theorem due to the sym- 
metry X°- — + X'^ + c"l, and therefore protected from 
quantum corrections 2]. However, the equation ([7]) for 
the embedding is modified at one loop. 

Let us focus on the first term J ct^xVC , which is 
essentially the vacuum energy due to zero-point fluctu- 
ations. In GR, it amounts to a huge contribution ~ L^ 
to the cosmological constant, which must be canceled 
by an extremely fine-tuned bare cosmological constant 
in order to reproduce the small value 0{meV) in the 
LCDM model. This is the well-known cosmological con- 
stant problem, which persists even in models with TeV 
scale supersymnictry. 

We claim that this problem is resolved here. To see 
this, note that 

\Gp.{x)\ = \g,U^)\ (14) 

independent of 9P'^{x). Hence the variation of the vac- 
uum energy term 

S J dt^xVC - j d'^x^g^'^'Sg^.u - J d^x^^^ScI)' Ag(j)^^j 
vanishes for harmonic embeddings 

Ag(j)' = 0. (15) 

Therefore the term J d^xVTj A"* has a different physical 
meaning here: It should not be interpreted as a cosmo- 
logical constant, but as a brane tension. In particular for 
harmonically embedded branes (i.e. minimal surfaces), 
the precise value of its coefficient ~ A^ is irrelevant and 
drops out from the equations of motion. For example, 
flat Moyal-Weyl space is a solution even at one loop, 
without fine-tuning A. Thus minimally or harmonically 
embedded branes (w.r.t. g^ii>) are protected from the cos- 
mological constant problem [3]; they are in fact stabilized 
through the vacuum energy. The crucial difference to 
general relativity is the parametrization of the geometry 
in terms of "tangential" 6'^" and "transversal" rather 
than a fundamental metric. Gravitons originate from 
fluctuations of 9^^'^ and are also blind to / d'^xVC . 

Combining ([7]) and , we see that harmonic embed- 
dings with gpi, ~ G^jy solve the vacuum e.o.m. of the 
matrix- model including quantum corrections (jl3[) . The 
presence of matter will lead to deviations from harmonic 
embedding, however these corrections are suppressed by 
factors 0{-^) where p is the energy density. Thus 
Ag4'^ = is valid as long as the brane-tension ~ A^ 
dominates the energy density of matter and curvature. 

Robertson- Walker geometries: Now consider cos- 
mological solutions of FRW type, and assume that Ai — 
0{TeV) to be specific. The Friedmann equations are 
replaced essentially by the requirement of minimal em- 
bedding with gpi, = G/^ii^. This should be valid to 



3 




FIG. 1: Embedding of universe in , for b = 1, m = 5 



high accuracy except possibly in the very early universe 
when the energy density becomes of order O(TeV). The 
FRW geometry has the general form 



ds^ = -dt^ + a{tfdY? 



(16) 



where 



1 — kr"^ 



dr^ + r^dQ^ = dx^ + Sixfdn^ (17) 



is the 3-dimensional metric with uniform curvature; here 
S{x) ^ r = (sinx, X, sinhx) for k = (1,0,-1) respec- 
tively. The following provides harmonic embeddings of 
FRW geometries for k = ±1 in K^*^, generalizing Ji]: 



where 



/ / S{x) sin6'cos(/J 

^ ^ 5(x)cos6' 

\C{x) 


V Xc{t) 



(18) 



and C(x) = (cosx,coshx) for k = (1,-1) respectively. 
Here r]ab in © is riab = diag(+, ±, -) for fc = 1, or 
Vab = diag(+, ...+,-, +, +) for k = -1. Then 

gf,^ dx^^dx" = r]ab dx'^dx'' ^ -c{t)dt^ + a{tfdT.'^, 

c{t) = k{xc^ -a^i;^ -d^). (19) 



This can be brought into FRW form (fTB]) using a change 
of variables ^ = \/'c(Tj. Equivalently, we can choose t 
to be the "proper" time variable r, such that a{t) is the 
usual FRW scale parameter; then c(i) = 1, and 



2„;,2 



k. 



(20) 



Note that the matrix coordinates and in particular Xc 
have no physical meaning from the brane point of view; 
they are determined by the requirement of harmonic em- 
bedding Agx" = 0. Due to the symmetry, it is enough 
to show that 



Ag{TZ{t)S{x)cos9) 



= AgXc- 



This leads to 



3- (d'^ + k) + d- ■ip'^a 
a 







5-00 -I- ipa — 

Xr = 0. 



„1 . . 

3— ax^ 
a 



(21) 

(22) 
(23) 
(24) 



These equations can be integrated as follows: 

(d^ + k)a^ + b^a^^ = m — const 

ip = ha^^ , h — const > 
a^Xc ^ d = const; (25) 

the last equation is in fact a consequence of ([20]) . which 
gives d = \frn and hence m > 0. This leads to 



i/2 = ^ = -fo^a-i" 



-8 



- K- (26) 



~3ma 



-10 



(27) 



For small a and 6 7^ 0, we find a ^ which signals 
inflation. For large a, it follows that a < 0, d ^ 0. The 
case 6 = was obtained before in [9|] . 

For fc = 0, one finds an unrealistic age of the universe 
of order 4.5 • 10^ years (assuming small 6), hence we will 
not pursue this case any further. It turns out that fc = — 1 
is the (near-) realistic case as discussed below. 

To complete the solution of emergent gravity we need 
to find a Poisson structure Q^^ which satisfies PT|) such 
that G 

— 9[i,v Clearly there exists no homogeneous 
and isotropic non-degenerate 0^'^ . However, e~'^ should 
at least be spatially homogeneous, because it determines 
the nonabehan gauge couphng Consider fc = — 1. 



dt 



We introduce t(t) through ^ 
metric in the form 



dt 



and write the FRW 



dsl^%-{^dT'' +dr' ^r'd??) (28) 
t^ 

where t — tcosh(x), t = tsinh(x). In particular, for 
a(t) = t we recover the well-known fact that the Milne 
universe is flat. Then the (complexifled) symplectic form 

6»-i = O-^dx'' A dx" = idr A dxi + dx^ A dx^ (29) 

in these flat coordinates is closed and («-) self-dual, hence 
(fTTj) is satisfied; recall that this equation applies even at 
the quantum level. Moreover we obtain G^y — 9^i^i 

and 



a" 



I 



(30) 
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In particular, e~'^ ^ 1 for large t as the geometry ap- 
proaches that of a Milne universe. Hence the gauge cou- 
pling is approximately constant for large t, but has a 
non-trivial time evolution in the early universe. 

Emergent cosmology: The cases fc = 0,+l imply 
a too short age of the universe given the present Hub- 
ble parameter. Therefore we focus on the case k — —1. 
Then d ^ 1 for large a, and the deceleration parameter 
9 = — f§ ^ 0. For large t, the time evolution approaches 
that of a Milne universe, which is in remarkably good 
agreement with observation [llj . The age of the universe 
is found to be -g- « 13.9 • 10^ years, the time evolu- 
tion of a{t) in the ACDM model at present being tangent 
with the evolution in the Milne universe Moreover, 
the main observational constraints including the acoustic 
peak in the CMB background and the type la supernovae 
data appear to be consistent with an interpretation in 
terms of a Milne Universe [ll|. While this geometry is 
excluded within GR, it makes perfect sense within emer- 
gent NC gravity. In view of the unreasonable fine-tunings 
in the presently favored ACDM model, this certainly de- 
serves a more detailed investigation. 

Inflation and big bounce: The scaling parameter is 
determined by 

d = \/-b^a-^ +ma-^ + 1. (31) 

For 6 7^ 0, denote with oq the (positive) root of the ar- 
gument, which is the minimal "size" of the universe. We 
fix the origin of time by a(0) = ao, and define ti by 

d = ^ a{ti) = \—=ai (32) 
V om 

using ([27| . Expanding l + ma~^ — b^a^^ — p{a — aa) + ... 

around oq where p = — \a=ao shows an 

inflation-like phase 

d ~ \/p(a^^ao), a{t) -|- oq, 

which ends at a{ti) = where d = ("graceful 

exit"). A typical evolution of a(t) and the corresponding 
Hubble parameter is shown in figure [2] Assuming < 
6'^, we have approximately 

„ ,1/4 pb^^ 
ao V 3 y/m 



Note that the requirements for "successful inflation" such 
as a large number of e-foldings will be greatly relaxed 
here compared with standard cosmology; this should be 
addressed elsewhere. For 6 > 0, it is obvious that the 
time evolution should in fact not start at t = 0, but be 
completed symmetrically as 

a{-t) = a{t), ij{-t) = -ij{t), x,{~t) = -x,{t) 

corresponding to a "big bounce" rather than a big bang. 




FIG. 2: Evolution of a{t) and H{t) for m = 5, 6 = 1. 



We conclude that Yang-Mills matrix models admit cos- 
mological solutions which are in remarkably good agree- 
ment with observation. The type la supernovae data are 
accommodated without any flne-tuning and without in- 
troducing dark energy. The solutions should be valid as 
long as the quantum-mechanical vacuum energy domi- 
nates the energy density due to matter or radiation, i.e. 
up to epochs with TeV-range temperature in the case 
of TeV-scale supersymmetry. Some modifications could 
arise from compactification in extra dimensions (leading 
to interesting low-energy gauge groups, effective scalar 
fields & potentials as in [141), or soft SUSY breaking 
terms (e.g. a mass term) in the matrix model. However 
it is unlikely that this would change our main conclusion, 
which is a Milne-like evolution after a big bounce and an 
inflation-like phase. If confirmed this would resolve the 
cosmological constant problem. 
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